Quantum Brownian motion in ratchet potentials: duality relation and its consequences 
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Quantum Brownian motion in ratchet potentials is investigated by means of an approach based 
on a duality relation. This relation links the long-time dynamics in a tilted ratchet potential in the 
presence of dissipation with the one in a driven dissipative tight-binding model. The application to 
quantum ratchets yields a simple expression for the ratchet current in terms of the transition rates 
| in the tight-binding system. 
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Brownian motion in ratchet potentials Q has attracted a lot of interest. One reason is the fact that ratchet 
systems, i.e. periodic structures with broken spatial symmetry, present the property of allowing transport under the 
00 ■ influence of unbiased forces. The interest has recently grown with the transfer of the problem in the quantum regime. 
7-H \ There, the description of dissipative tunneling [J] presents a theoretical challenge which was tackled in relatively few 
works [E 0, IE IE IE 0], while first experimental realizations were reported [E El • After the semiclassical work Q , 
further progress towards the theoretical description of quantum ratchets involved modeling in terms of tight-binding 
systems 0, Q or in terms of a molecular wire Other available methods include perturbation theory or a 
quantum Smoluchowski equation Q . Most of these methods 0, S IE 01 are restricted to the regime of moderate-to- 
strong friction. 

The approach discussed in this article originates from works Il2t IT^l Il4. Ilfil| on quantum Brownian motion in 
a tilted sinusoidal potential, which led to a duality relation for the mobility of the system considered with the one 
of a driven dissipative tight-binding model. This idea emerged first in Ref. |ll| , where the linear dc mobility at zero 
temperature was considered. In that work the interest was focused on the occurrence of a diffusion-to-localization 
transition in the system with increasing dissipation. These results were corroborated by means of renormalization- 
group methods |l2j |. The duality relation was extended to the nonlinear dc mobility at finite temperatures in Ref. |13| . 
where an extensive physical discussion as well as important milestones of the proof were given. It was subsequently 
3 ■ applied to the investigation of the current- voltage characteristic of small Josephson junctions |14| . Later, an identical 
duality relation for the linear ac mobility was obtained by a different approach in the frame of linear response HE. In 
particular, that work went beyond the case of a strictly Ohmic dissipative bath considered in Ref. Il3l. and included the 
case of an Ohmic bath with finite cutoff frequency as well as sub and super-Ohmic baths. In Ref. |1(t| . the formalism 
of Ref. HE| was generalized to arbitrary ratchet potentials, i.e., periodic potentials of arbitrary shape. Moreover, the 
duality relation was extended to the average position of the quantum particle at long time. In the present article, we 
give detailed proofs and discussion of these last results. The duality relation is obtained in its most general form, in 
terms of the generating function out of which the average position, the mobility, and other dynamical quantities can 
| be extracted. Its domain of validity includes weak dissipation and nonlinear driving. However, our demonstration 
remains restricted to the case of a strictly Ohmic bath. 

The article is structured as follows: The problem is defined in Section The duality relations and their discussion 
are presented in Section [HI The application to the evaluation of the ratchet current is developed in Section llTTI The 
details of the proofs are found in Section llVl 
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I. STATEMENT OF THE PROBLEM 

We consider the Hamiltonian Hr of a quantum particle of mass M in a one-dimensional periodic potential V(q) 

O ■ tilted by a force F, 
O 
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The potential can be any function of periodicity L, and is fully characterized by the amplitudes V\ and phases tpi of 
its harmonics in the Fourier representation 



V(q)=^V l cos{2-Klq/L-y l ). (2) 
1=1 

In order to investigate quantum Brownian motion, we have to let the particle interact with a dissipative thermal 
environment. This is modeled by the standard Hamiltonian iJe of a bath of harmonic oscillators whose coordinates 
are bilinearly coupled to the system coordinate q 
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The bath is fully characterized by its spectral density J{lo) — (tt/2) ^a=i( c a/' m a LL >a)8(uj — Lo a ), defined in terms of 
the masses m a , frequencies u> a , and coupling strengths c a of the oscillators. We consider an Ohmic spectral density, 
i.e. linear J(w) ~ r\ui at low frequency u>. The viscosity coefficient r), together with the particle mass M , defines the 
time scale of dissipation 7" 1 = (rj/M)^ 1 . 

The information on the system dynamics is contained in the reduced density matrix p(t) = Tre W(t), obtained 
from the density matrix W(t) of the system-plus-bath H = Hh + Hr by performing the trace over the bath degrees of 
freedom. The diagonal elements P(q, t) = (q\p(t)\q) of the reduced density matrix suffice in order to evaluate, e.g., the 
evolution of the average position (q(t)) — Tij\{qp(t)}. It turns out to be very powerful to work with the generating 
function defined as 

P(A,t)= / dqe^P(q,t). (4) 



The normalization of the reduced density matrix implies the property P(X — 0,i) = 1. The name and utility of 
the generating function come from the fact that its derivatives generate expectation values of powers of the position 
operator 
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In particular, we would like to investigate quantum Brownian motion in a ratchet system. Apart from special 
configurations {Vi sin (cpi — l(pi) = V7} of the amplitudes Vi and phases (fi, the potential J5J is spatially asymmetric 
and describes such a ratchet system. The ratchet effect is characterized by a nonvanishing average stationary particle 
current wg 2 = limt_, 00 t^ 1 J Q dt'v(t') in the presence of unbiased time-dependent driving F(t'), characterized by 

limt-xxj t^ 1 J Q dt'F(t') — 0, and switched on at time t' = to- In this article, we report a method to evaluate the 
stationary velocity v^ c (F), which is simply obtained by time-differentiation of (q(t)), in the biased situation (JJJ of 
time-independent driving F. The ratchet current in the presence of unbiased bistable driving switching adiabatically 
between the values ±F is obtained through the relation = v™ c (F) + v^ c (—F). 

The diagonal elements P(q, t) of the reduced density matrix can be obtained by real-time path integrals techniques 0, 
IT^ I . At initial time t' — to, we assume a preparation in a product form where the bath is in thermal equilibrium with 

the system W(to) — p(to)e~ f3 ^ ' [Tre e _/3 ^B 'j -1 . In H B °\ the system operator q is replaced by an initial position go- 
The bath temperature is fixed by T = l//3fcB- This leads to the expression 



P(q,t) = (q\ P (t)\q) = J d Ql J d^{ qi \p(to)\^)G(q t q, q U t) (6) 
with the propagating function 

G(q f ,q' f , qi ,qi,t)= Vq ( * V*q'A[q, t]A* [q', t]Fpy[q, q', t] (7) 

J qi J q'. 

given as a double path integral on the paths q(t') and q'(t'). For the continuous coordinate q and a Hamiltonian of 
the form (J), the path integral stands for [17J, ch. 2] 
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where the time interval t — to has been sliced in Ni intervals of length At = (t — to)/Ni and the path q(t') has been 
discretized into the set of values qk = q{to + fcAr) for k = 1, . . . , Ni — 1. The boundaries of the path integral remind 
of q(to) = qi and q(t) = qj. The propagator A[g,t] reads 



A[q,t] = exp {- [ dt' 
i h Jt n 



M 

-q(t'Y-V{q(t')) + Fq{t>) 



(9) 



The Feynman- Vernon influence functional may conveniently be written as Fpy[q,q',t] = exp{$pvK, X, t]}, in terms 
of the phase 
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|Mj(0) / dt'£(t')x(*') + T«(*o) / dt'CCtOM^t'-to), (10) 
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which depends on the difference £(t') = g(t') — g'(t') and average x(t') = [9(t') + 9'(^')]/2 of the two paths. It 
induces nonlocal-in-time Gaussian correlations between the paths. The functions Lr(t) and L\{t) denote the real 
and imaginary part of the bath correlation function 



1 f°° 
L{t) = - / dujJ(uj) 
Jo 



coth 



cos(wr) — isin(wr) 



(11) 



Integrating the imaginary part yields the function M\(t) — duij(uj) cos(uit)/ttu>. 

Our goal is to evaluate the generating function JIJ. In general, the nonlinearity of the potential V(q) prevents the 
evaluation of the path integrals. However, an exact expansion of the contribution of the potential in the propagator A\q] 
makes it feasible. In the following, we shall present this method, which was first introduced by Fisher and Zwerger |l3| 
for a sinusoidal potential. 



II. RESULTS AND DISCUSSION 

In order to expand the contribution of the potential, we introduce a "charge" a taking values in the set {±1, ±2, . . .}, 
and corresponding amplitudes defined as 



A a = ye ivv for g > 0, 

This allows to rewrite the potential (01 as a simple sum of exponentials 

V(q) = &*c~ 2 ™ q/L - 
<r=±l,±2,... 



a;. 



(12) 



(13) 



Using this expansion, one can demonstrate that (see Section TlV A|l 
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where we have introduced n charges Oj and corresponding times tj, as well as the function p(t') = 
(2nh/L) Y^j=i a j&{t' — This expression can be substituted in the propagator A[q, t]. We have to do the same 
for the second propagator A*[q',t], using the complex conjugate of Eq. I|14fl with a new set of n' charges a'j, and 
corresponding times tj,, and a new function p'(t'). The product A[q, t]A*[q', t] may then be conveniently rewritten in 
terms of the difference £(t'), respectively average path x{t')i 



A[q,t}A*[q',t] = 
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dt' 
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(15) 



The gain of this expansion is that the paths now enter at most quadratically in the argument of the exponential. 
Eventually, the path integrals will become Gaussian integrals. The price paid is the emergence of a series expression. 
We have introduced a compact notation for the sums, products and integrals involved 
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(16) 

Performing the now Gaussian path integrals, and after a long calculation described in Sections II V Bl ITVEl we obtain 
our main result 



P{X,t) 
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(17) 



As indicated by the relation symbol, this result is valid within some approximations, as shown in the full derivation of 
Section llVl The validity regime of (|T7jl is discussed in detail below Eq. The functions y s h(t') and x s h(t') denote 

the difference y s h(t') = q s h(t') — <?sh(*')> respectively average x s h(t') = [q s h(t') + q' sh (t')]/2, of the step-like paths 



q sh {t') =LJ2<T j [6{t' - tj) - 1] + q sh (t) 

3 = 1 

n' 

3' = 1 



(18a) 
(18b) 



The step heights are multiples of L = 2Trh/r]L. The paths end up at the same value q s h(t) at the final time t. Therefore 
the difference path ends up at y s h(t) = 0. The quantities 



AC = y S h(t) - y S h(to) = L 



Ax = x sh (t) - x sh (io) = 
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specify the initial boundary condition for the two paths. The primed sum 
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(19a) 
(19b) 
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where 6 denotes the Kronecker symbol, is thus restricted to the configurations for which the difference path starts 
at y s h(to) = 0. The influence phase <I>py [y s h , x s h , t] is defined as in (|10fl . provided that the spectral density J(ui) 
entering the correlation functions L(t) and M\(t) is replaced by the new spectral density 



Jtb(w) 



J{u) 



1 + ( W / 7 ) 2 



(21) 



Finally, the auxiliary function N{t) is discussed in Section llV El 

The justification of these notations appears when one considers the generating function Ptb(A, t) of a driven tight- 
binding model given by the Hamiltonian 
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(22) 



m—l l- 



The couplings A rn are precisely the one introduced in l|12(l and involved in the boxed sum. The spatial periodicity 
of this tight-binding model, which can be read in the position operator (Jtb = -^X^-oo is precisely the height 

unit L of the steps of the paths Ijl8(l . This tight-binding model is bilinearly coupled to a different bath of harmonic 
oscillators characterized by the spectral density Jtb(^) given in Q2ip. This spectral density is still Ohmic, with the 



same viscosity coefficient 77, but now presents a Drude cutoff at the frequency 7 set by the dissipation of the original 
model. The system is initially prepared in the state /5tb(^o) = |'o)(fo| with IqL = q s b(to)- hi this situation, the 
generating function reads 



Ptb(A^) = 



(23) 



which bears a clear structural resemblance with (|17l) . The A 2 -terms are absent of (|23|) . but they do not play any role 
as far as one is interested in the average position (q(t)} [see Eq. JJJ]. One also notices that A\ comes with an opposite 
sign in the two expressions. 

The link between the original model Q &ud the tight-binding model ()22)1 can be pushed further. The normalization 
of the generating function P(X = 0, t) = 1 yields the identity 



iF 



exp j$^[y sh ,x sh ,t] + y y dt'y sh (i'; 
starting either from (| 1 T|l or 1)23)1 . Differentiating with respect to i 7, yields the set of non-trivial identities 



(24) 



= 



£' 



dt'y ah (t') 



exp j$^[y sh ,x sh ,t] + y jf dt'j/ sh (t')|, 



(25) 



for any fc = 1, 2, . . . , 00. 

Another important result can be obtained from the relation (|17|) by evaluating the average position (q(t)) 



[dP(X,t)/dX] 



Using the identities (|24|l and 1)25)1 with k = 1, one gets 



(«(*)) ~ (<f(M) 



(p(t )) , F(t-to) 
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(9TB (i)) 



TB, 



(26) 



where (q(to)) = Tr^ {<7p(io)} and (p(io)) = Tr^ {p/3(to)} denote the position and momentum of the initial preparation 
of the ratchet system. The last term of 1)26)1 is the average of the position operator qtb in the driven dissipative 
tight-binding model ffify. initially prepared in the state /3tb(*o) = |0)(0|. It can be obtained from Eq. with 

(<Ztb (t)) — [9Ptb (A, t)/dX] and Iq — 0. It comes with a minus sign due to the minus sign in front of Ax in 1)1711 . 

The duality relation for the position (|26() is a very useful result for quantum ratchet systems. We shall discuss this 
application in Section IlIII 

The duality relations l|17l) and (|26|) are approximate results, as denoted by the relation symbol. As derived in this 
work, they arc valid when the following conditions are simultaneously met: 

i) Long-time dynamics: The measurement time t — to should be much longer than the time scale I/7 set by 
dissipation. This can be easily controlled experimentally. 

ii) Rare transitions limit: The terms e -7 ^' 1 -* ), e ~ 7( * -ttl) , e -" B (* tr ~* ), and e -" B (*~ ttr ), with uj b = 2nk B T/h, should 
be negligible with respect to 1, when t tr equals any of the times tj,t'j,. These times, which are integration 
variables involved in the boxed sum, are the transition times in the double path integral representation 1)23(1 of 
the generating function of the tight-binding model. Therefore, this approximation corresponds to neglect, in 

g 1 ) after 

t. It will therefore be valid when the transitions in the 



the boxed sum, the contributions from the paths which involve transitions on a time scale max(7" 
the initial time t' = to or before the final time t' ■■ 
tight-binding model are rare on a time scale max(7~ 
the temperature. 



This condition is controlled by the dissipation and 



Furthermore, in our derivation we have used a strictly Ohmic spectral density J(u>) — rjuj. In this case, the 
function Mi(r) takes the simple form M\{t) = r)5(T), and the divergence of L™(0) allows to restrict the configuration 
sum to its primed version. A physically more realistic situation would be to consider the Ohmic spectral density J{uj) — 



r/uje 



with finite cutoff frequency uj c 



We do not know to which extent these restrictions are specific to the method that we have used in order to derive 
the duality relations. An equivalent duality relation for the mobilities [see Eq. ()27)l ] has been derived in Ref. ^5| in 
the frame of linear response for a sinusoidal potential. It is interesting to notice that the derivation presented in that 
work does not require the restrictions to the rare transitions limit and to a strictly Ohmic bath. A more generalized 
version of l|2*7| has also been obtained for a much broader class of spectral densities, including sub-Ohmic and super- 
Ohmic ones, and in the case of time-dependent driving. We currently do not see any problem of principle in order to 



generalize our demonstration for a general form of the spectral density, and this is the subject of work in progress. 
However, we do not see how to avoid the restriction to the rare transitions limit in our derivation. Furthermore, we 
do not know how to generalize the identities 1)25(1 . which we have used in our proof, for the case of time-dependent 
driving. These remain open questions. 

Let us now give some interpretation of the results. The duality relation (|17|l for the generating function is not very 
useful in itself, but very powerful in order to generate useful results. It links the dynamics of the two systems Q 
and (J22J- The precise relation between the two systems is specified by: i) The relation (|12|) between the harmonics 
of the potential of the original system and the couplings in the tight-binding system; ii) The relation L = 2Trh/r]L 
between the spatial periodicities L of the original system, and L of the tight-binding system; iii) The relation l|21() 
between the spectral densities of the baths of harmonics oscillators coupled to each of the two systems. In the original 
system, the relevant dynamical parameters are captured by the dissipation parameter a — r\L 2 /2ttH and the energy 
drop per periodicity length e = FL. Due to the change of periodicity length, these parameters become a = 1/a 
and e = e/a in the tight-binding system. Thus, weak dissipation in one system maps to strong dissipation in the 
other one, although the viscosity rj in the spectral density does not change. 

The duality relation for the average position l(26|l is an example of a useful result which can be extracted from l(17|l . 
There, the relation between the average positions (q(t)) and (qtb (^))tb in the two systems, which holds at long time, 
is explicit. The asymptotic dynamics covered by this result is usually described in terms of the nonlinear mobility 
/.t = limi_ >00 (g(t))/i ? . Accordingly, the duality relation (|26|l may be rewritten in the form 

H{a,e) ~ /i - HT B {l/a,e/a), (27) 

where //o = 1/?/ is the mobility of the free system, V(q) = 0. In the special case of a sinusoidal potential, this relation 
was already obtained in for the dc mobility. As mentioned above, it has also been derived in 15] for the linear 
ac mobility in a sinusoidal potential. 

The second derivative of the generating function P(X,t) with respect to its parameter A yields the variance (q 2 (t)), 
which gives information about diffusion and current noise. It would thus be natural to try to extract from 117JI a 
duality relation for this quantity. However, the result diverges, because the quantity N(to) involved in the right-hand 
side diverges for the strictly Ohmic spectral density J(w) = rjuu considered in the derivation. In order to get results 
on diffusion and current noise, we shall thus have to go beyond this approximation and allow for an Ohmic bath 
with finite cutoff frequency. As mentioned above, we do not see any problem of principle in order to generalize our 
demonstration to this situation. 



III. APPLICATION: EVALUATION OF THE RATCHET CURRENT 

In this section, we shall discuss the application of (|26|1 to evaluate the current in ratchet systems. By time- 
differentiation of (q(t)), given on the left-hand side of (|26J) . one obtains the stationary velocity v B ? c (F) in the biased 
situation of time-independent driving F. As discussed above, the ratchet current in the presence of unbiased bistable 
driving switching adiabatically between the values ±F is obtained through the relation v^ 1 = v^ c (F) + v^ c (—F). 
Our task is therefore to evaluate the right-hand side of l|26(l , in particular the average of the position operator qtb in 
the driven dissipative tight-binding model l|22(l . initially prepared in the state ptb(^o) = |0)(0|. It can be obtained 

from the generating function Ptb(A,£) of the tight-binding system through (<7tb(0) = [dPpB (A, t)/dX] 

x=o 

The generating function can be obtained from with Iq = 0. However, this formula is not the most suitable 
in order to get the long-time behavior required in i|26|) . For this purpose, we go back to the population P ; TB (t) of 
the site \l) of the tight-binding system at time t, after a preparation in the site |/o) = |0) at initial time to. These 
populations are the diagonal elements of the reduced density matrix of the system ()22|l coupled to its bath of harmonic 
oscillators of spectral density Jtb(w). They are related to the generating function by Ptb(A, t) = X^-oo e AiL P; TB (i). 
The real-time path integrals techniques which led to (|23(l yield for the populations an expression which can be rewritten 
in the form of an exact generalized master equation 



r-t °o 

i?»(t)= / df £ K^{t-^P^{t'). (28) 

Jt ° l' = -oo 

For the rather technical discussion of the expressions for the kernels if m (r), we refer to Refs. 0,0]. The generalized 
master equation can be easily rewritten in terms of the generating function, yielding 

P TB (X,t)= [ &t'K{\,t-t')P TB {\,t'), (29) 

J to 



with K(X, t) = X)m=-oo ^ XmL K m (r). As shown in Section TlV Fl the solution of this equation behaves at long times 



as 



Ptb (A, t) ~ expJ(t-to) Vr m (« 



e Ami - lj | . (30) 

The transition rate T m from a site |Z) to a site \l + m) is obtained from the corresponding kernel by T m = J °° dri^ m (r). 
From this result we derive easily 

(£rn(t)) - (t-t )iymr m . (31) 

Plugging this in the duality relation (|26[1 and differentiating, we obtain the stationary velocity v^ c (F), respec- 
tively v^ c (—F), in terms of the transition rates T m (F), respectively T m (—F), 



v 



DC (±F) = ±--LJ2 mr m (±F). (32) 



The ratchet current in the presence of adiabatic bistable driving reads accordingly 

v r =- L J2 m i r ™( F ) + r ™(- F )i • ( 33 ) 

This result shows that the ratchet current in a system characterized by the potential (J2J) is very simply related to 
the transition rates in the tight-binding model (|22[l . As the duality relation from which it is derived, it is valid in 
the rare transitions limit T m (±F) <C min(7, 2 / nk^T/K) and for a strictly Ohmic bath characterized by the spectral 
density J(lu) — rjuj. 

The ultimate task is thus to evaluate the tight-binding transition rates T m (±F). The expression obtained by 
real-time path integrals techniques consists of numerous contributions which can be classified with respect to the 
number N of transitions in the double tight-binding path which they involve. We denote by T^m\±F) the sum 
of all contributions involving N transitions and call it the Ath order transition rate. The total rate follows from 
r m (±F) = X)jv=2 \±F). The derivation of the explicit expressions are discussed in Section HV Gl At second 
order, we find, for m/0, 

r (2)(F) = ^ ATe -~raHL-/H)Q{r) + - ml (FL/H)r^ ( g 4) 



The twice-integrated bath correlation function is defined as 



s 1 t°° , Jtb{u) \ , fhw/3\ r . .. . . . . 

Q(T) = — / dui - — coth 1 — cos(wr) + lsm(cjr) 

k Jo v 2 I V 2 / 



(35) 



with the spectral density Jtb(w) given in 121|) . We remind that the couplings A m are related to the harmonics of 
the potential (0) through Eq. i|12|) . One sees that the phases <p m of the potential, which are identified to the phases 
of the couplings A m , do not come into play at second order. 

At third order, there are contributions which involve only A±i and A±2- They read 

iffi [112](F) = 2 * Al * ^ Im ( / dTG[ 3) | [112](r)e im ^ T - isgn(m)vil2 l , (36) 



fr 3 

for m = ±1, ±2, with the functions 

/"°° - 2 r - 2 ~2 n 

Gl 3) [112](r) = - / dpe-^^rt e ^[-2Q(r+p) + Q(r+2p)] +e i-[-2Q(r-p)+Q(r-2p)] ^ 

Jo L 

r°° -2 

g( 3) [112](tH / dpe^[«^- 2 ^ T+ ^)- 2 ^-^)]. (37b) 
Jo 

The phases ip m of the potential enter through the unique phase difference P112 — f2 — 2ip±. 
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FIG. 1: Stationary velocity (a) and ratchet current (b) as a function of temperature, for different values of the driving amplitude. 
Weak dissipation is chosen with a = 0.2. For all curves, the point at the lowest temperature is evaluated for ksT/AV — 0.1. 
The lines are guides for the eyes. 



There are also contributions involving A±i, A±2> an d A±3. They may be rewritten 



r^)[123](F)= 2|All L A2||A3 'Im 



drGf2[123](r)e im -TT r - isgn ( m ^ 123 



for m = ±1, ±2, ±3, with the functions 



g[ 3) [123](t) 
G^ 3) [123](t) 



Ape 



K-p) 



V[-3Q(T+2p)+2Q(r+3p)] 



+ e n 



-3Q(T-2p)+2Q(i--3p)] 



(38) 

(39a) 
(39b) 
(39c) 

and the phase difference tp\23 = (fia — <f2 — <P\ = (</>3 — 3y>i) — {f2 — 2</?i). 

If the potential (J2J) sustains at most the first three harmonics, meaning that the couplings A m are all for \m\ > 3, 
there are no other contributions to the third-order rates, which are then given by 



dpe 



-^Q(-p) 



4r[-6Q(r+ip)+2Q(r+fp)] + ^ [-6Q(r- ip)+2Q(r- § p)] 



G^ 3) [123](t) = / dpeTQ(")[ e x[-6Q( T 4f)-30(r+Ip)] + e ^[-6Q<V+ip)-3Q(T-§p)] 



r^( J F) = rW[ii2](F) + rS ) [i23]( J F). 



(40) 



A very interesting feature of the expressions for the second-order and third-order rates is their explicit and simple 
dependence on the moduli and phases of the couplings A m , related to the amplitudes and phases of the original poten- 
tial (0). In particular, this dependence has direct consequences for the ratchet current (|33|l . Firstly, one can explicitly 

(2) (2) 

verify that the ratchet current vanishes for symmetric potentials. The second-order rates satisfy Tfn (-F) = r_} n (F), 
therefore they do not contribute to the ratchet current. Likewise, one can see that the third-order contributions to 
the ratchet current, which are the dominant ones, are proportional to sin(<^>2 — 2yi) or sin(<p3 — 3(pi). Therefore, 
they vanish for spatially symmetric potentials, characterized by {V/sin(<p; — Itpi) — V7}. Secondly, this simple 
dependence on the potential parameters should be accessible in experimental realizations where the potential can be 
tailored, as, e.g., in arrays of Josephson junctions |Tfij . 

We see that the complexity of the transition rates increases with the order N. In tight-binding models with large 
dissipation parameter r]L 2 /2irh and/or high temperature, neglecting higher orders is known to be a good approxima- 
tion j2j. 

Finally, we discuss the behavior of the stationary velocity l|32|l and ratchet current l|33|l as a function of the driving 
amplitude, the temperature, and the dissipation strength for a potential sustaining two harmonics only. The rates are 
evaluated up to third order, using the expressions and iJSBJ. The outcome is shown in Figs-HandE] With V\ = 4V 2 , 
the untilted potential, obtained from (J2J, has a barrier height AV^ = 2.2V\. Another energy scale relevant for the 
dynamics is the confinement energy E c { — 2ttH 2 /ML 2 of a particle of mass M confined within a length L. The 
dynamical regime can be specified by the ratio between these two quantities. With AV/E C { = 0.264, the typical 
action of the particle, S = V2MAVL 2 = ^/AV/E c f7r x h, is about Q.3h. The dissipation strength a — r/L 2 /2ttH is 
the ratio between the energy scale fvy associated with the dissipation rate 7 = T)/M, and the confinement energy E c f. 
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FIG. 2: Stationary velocity (a) and ratchet current (b) as a function of dissipation strength, for different values of the tempera- 
ture. Weak driving amplitude is chosen with FL/AV = 0.01. In this regime, the stationary velocity (a) is a good approximation 
of the linear mobility, given in units of no = 1/rj. The lines are guides for the eyes. 



Another measure of the dissipation strength is given by the ratio between 7 and the classical oscillation frequency Qq = 
2w y/Vi / M L 2 in the untilted potential. With the value a = 0.2 used in Fig. [IJ this ratio, j/Clo = a / 2ttV\ / E c f, is 
about one fourth, which corresponds to weak dissipation. 

The behavior of the stationary velocity as a function of the temperature for different values of the driving amplitude 
and weak dissipation a = 0.2 is shown in Fig.GJi. The first curve at low driving FL/AV = 0.01 (circles) is also a good 
approximation of the linear mobility, given in units of (j,q = 1/V The first value at low temperature kBT/AV = 0.1 
is very close to the velocity vq — F/rj obtained in the absence of the potential V = 0. Upon increasing temperature, 
the stationary velocity first decreases, before increasing again above a crossover temperature k#T* / AV » 2. This 
behavior is characteristic of weak dissipation. This crossover temperature is lower for the third curve, at larger 
driving FL/AV = 3 (up triangles). The ratchet current for the same parameter values is shown in Fig. ^d. A 
non-monotonic behavior with reversals is found. From these two figures, one also sees that the stationary velocity 
tends to vo for driving amplitudes or temperatures much higher than the potential barrier, and that the ratchet 
current vanishes correspondingly. A discussion of the behavior of the stationary velocity and of the ratchet current 
as a function of the driving amplitude for the same system can be found in Ref. [l(| . 

The behavior of the stationary velocity as a function of the dissipation strength for different values of the temperature 
and weak driving FL/AV — 0.01 is shown in Fig. E^. At weak dissipation, the stationary velocity (and equivalcntly 
the linear mobility) tends to the value vq obtained in the absence of the potential upon decreasing the temperature, as 
already seen in Fig.^,. This indicates delocalization of the ground state. At strong dissipation, one finds an opposite 
behavior. The few curves shown in the figure indicate that the transition between these different behaviors happens 
around a = 1 . The analysis of the power-law dependence of the transition rates at low temperature shows that this 
transition happens at a = 1 |l6[ . The ratchet current, shown in Fig. |2b for the same parameter values, explicitly 
reveals the occurrence of reversals as a function of the dissipation strength. 

In these numerical applications, none of the rates exceeds 0.047 and 0.05cjb, which means that the duality relation 
is valid for this system. However, one should keep in mind that only contributions to the tight-binding transition 
rates up to third order in tunneling have been evaluated. Higher-order contributions could become relevant, e.g., at 
strong dissipation, when the dual tight-binding system experiences weak dissipation a = 1/a. 



IV. PROOFS 



A. Expansion of the potential 



First, we shall demonstrate how the expansion (|13fl of the potential leads to the series expression (|14|) . Using the 
power series representation of the exponential function, we write 

exp <^ -- / At'V{q{t')) \=y2 dt n / dt„_! . . . / dti n 

I n J to ) n=0 ^*o J to J t j=1 

Due to the complete symmetry of the integrand in all the tj, the n integrals have been entangled, compensating the 
1 /n! factor coming from the series expansion. We now take advantage of the representation l|13|) of the potential as a 



-jV{q{t 3 )) 



(41) 



sum, introducing a charge (Xj and amplitudes A CTj . for each term V(q(tj)). By distributivity, product and sum can be 
exchanged, yielding a sum on configurations = Em=±i ±2 ■ ■ -So- n =±i ±2 ■ This leads to 



n 



En )«p^~ E^fe) 



(42) 



With the help of the function p(t') introduced below Eq. 1)14(1. this can be rewritten in terms of the path q(t') 



n 

3=1 



En(-%)-p{4//w)} 



(43) 



Substituting this expression in (|41(1 demonstrates l)14[l. 



B. Evaluation of the path integrals 

We now turn to the evaluation of the path integrals. We start from the expression for the propagating function • 
We rewrite the path integrals in terms of the difference respectively average path x(i'), with the boundary 

conditions & = g< — q[, Xi — {li + /% £/ = Qf ~ l'fi an d Xf — {if + Collecting the expressions l(TU}l and JTSJl, 

we have 



G(q f ,q' f ,q i ,ql,t)= £ jf ' ^ ^ ©* X exp{-,S R [£] - i5i[fl + y jf d^(i') 
+ i df X (f ) [-M|(t') + ^(f ) - [p(0 - + '-f [t(t) Xf ~ £(*o)x 



fi -/to Jto 

ft[e] = ^*d^(i')5[p(0 + />'(0]- 



HZ 

?! 



~£fXf + &Xi + 



with the definitions 



(44) 

(45a) 
(45b) 



For the imaginary part of the bath correlation function, we have made use of the simplification M\(t) = i]S(t) 
for a strictly Ohmic bath. We also have performed partial integration in order to remove the dependence on the 
derivative x(t') of the average path. 

We discretize the paths according to the procedure described below Eq. JHJ. The integrals in the argument of the 
exponential are discretized as in the following example 



to 



di'x(t') -Mi(t') + v £(t') - w) - P '{t')} = At ^ + nik - [Pk - p' k ] 



Ni-l 



k=0 



(46) 



and we evaluate the derivatives of the path £ with the difference formulae = (£,k+i — ^fe)/Ar and = (£fc+i 
+ ^fc_i)/Ar 2 . The terms involving Xk form the integral 



dx/c exp { -Xk 



-M'i k + rii k - [p k - P ' k ] 



2ttS [ - Xk 



-M£, k + r]ik - [Pk 



Pk\ 



(47) 



This 8- function allows to suppress the integral on The process is repeated for all values of k = 1, . . . ,Nj — 1. Some 
care has to be taken with the prefactor of in the argument of the <5-function, which will come as a denominator 
in front of the expression, and with the behavior of the path £(t') and its derivatives at the boundaries. One already 
sees that the result of the whole process, in the limit N\ — > 00, is that the path £(t') will be constrained to follow the 
solution y(t') of the differential equation 



-My(t')+r,y(t>) = p(t')-p'(t') 



(48) 



with boundary conditions 

2/(*o)=6, V(t) = €f- 

After having performed all integrals and taken the limit Nj — ► cx>, we obtain 

rt 



(49) 



G(qf,q' f ,qi,q'i,t) 



V 



2nh(l-e) 



J2 exp|-5 R [y] - iSi\y] + y j[ di'y(t') 



+ — fo(*)x/ - y(*oki] + ^ 



-Ox/ + &Xi + 2^ 



, (50) 



where we have introduced the notation e = e _7 (*~*°). The solution of the differential equation (|48|l with boundary 
conditions (|49|l can be written as y(t') = 2/hom(i') + J/part(*')> m terms of a solution of the associated homogeneous 
differential equation 



2/hom(0 = 



f 



1 - e -7 ^ - *')" 


+ i 4 " 


e 7(t'-*o) _ 




1 — e 





(51) 



and a particular solution of the differential equation 



3=1 



0(f - *i) 



l_ e 7(f-*i) 



[e - e~ 7( *J~ to:i ][l - e 7(Y_ * o) ] 



j'=i 



1 - e 



7(t'"t'/) 



1-e 



r £ _ e -7(*^-to)- 



[1-e 



7(*'-*o)l 



1-E 



(52) 



The periodicity length of the tight-binding model L = 2nh/rjL comes into play at this stage. 



C. Evaluation of the generating function 



The result (|50(1 for the propagating function can now be used to obtain the generating function. One has first to 
evaluate the integral J dqe Xq G(q, q, qi, q[,t) . Again, it will yield a 5-function for £j 



dq exp < q 



iM 



m + A 



with A£ defined in fljj) and A£ c = L \Y%=i <Tje~^< {ti ~ to) - £™,' =1 4e" 7(i >"' o) 
after transforming the variables qi and q[ of the initial integrals into and \ii 

ih{\ - e) 



(53) 



The generating function reads then, 



P(X,t) 



E 



d£i8 & - A£ c + A£ 



IF 



17 7 e 2 



A J exp j-S R [y] - iSrh] + — jf dttyf ) + ^£ 
x / dXi(Xi + e</2|p(to)|x< " 6/2) exp(-^-y(t )x, + ~&t< \- (54) 



In the position representation, an integral of the form J d\i(Xi\ ' \Xi) is a trace Tr R {-}. With the help of the eigenstates 
of the momentum operator, which are described by the wave-function (q\p) = e lqp ^ h /\ // 2nh, one can rewrite after some 
algebra 



dXiiXi + 6/2|p(to)|x< - 6/2) exp 



iM vtj I , 

-TVitolXi + -r&Xi t = Tr R i PW exp 



-iMy(io)+iij6 . , i& 



(55) 



In Eq. (|54|l . the integral on £j can be removed, provided that one substitutes everywhere in the integrand by the 
value given in the argument of the (5-function. In the path y(t'), this substitution yields y(t')\ £ i= _Af+A£ -ift(i-e)A/j7 = 
2/ sm ,A(i')> wi th 



ih\ r 

17 



1-e 



-7(*-*') 



■yem(t')- 



(56) 



FIG. 3: Typical smeared path j/ sm (i') (left) and sharp path y s h(t') (right), as defined in Eq. Ij57^ . respectively above Eq. 
The height of the steps are multiples of L and the edges of 2/ sm (t') are smeared on a scale I/7. 



The path y sm (t') — y S m.\=o(t') assumes the expression 

n 

y sm (t') =1^2*3 W - h) - 1] [l - e^'-^ 



iE4 [^'-^)-l] [l-e^) 



(57) 



The index sm stands for smeared, because y sm {t') is a step-like path whose edges are smeared on a scale I/7, as one 
can see on an example drawn in Fig. Collecting everything, one obtains the intermediate result 



P(X,t) = 



E 



Tr R < p(*o) exp 



-|A£ + Aj 9 + f-f (A£ - A£ c ) + (1 - e)Xj P - 



IF f vq 
x exp -j -5r[2/ 8TO ,a] - iSi[y BBl ,A] + J dty sm ,\(t) + — 



-A£ + A£ e -^-^A 



(58) 



D. Identification with a tight-binding expression 

The following task is to rewrite the generating function (|58|l in terms of the sharp path y s h(t') defined above Eq. Q18|l 
and represented in Fig. |3J The easiest way to understand the mechanism of this transformation is to work in Fourier 
representation. Technically, because these paths are defined in the time interval [to,t] only, one has first to continue 
them in the whole time axis by defining y(t') = \6{t' — to) — 6(t' — t)]y(t'), in order to be able to use the usual 
Fourier transform y(uj) — dt'y(t')e^ luJ ^ and take full advantage of the usual differentiation and convolution 
properties. One can then demonstrate the relation 



uj + vy 



AC - A^e 

7 



with the definition 



y s h,x( t ') = - 



9(t' - t ) - 6{t' - t) - S(t' - t ) 



1 - £ 



(59) 



(60) 



This can be done by considering the differential equation 1|48|) in Fourier representation. This relation means that, up 
to a boundary term, the Fourier transform of the smeared and sharp paths are related by a factor vy/(u> + vy). Let us 
also write the real part of the influence phase Srsjysm] in Fourier representation 



1 f°° 

SR[y sm \ = / du)J(u)coth(fw/3/2)y sm (u)y sm (-uj). 
J 



(61) 



It can be rewritten in terms of ysh(u>) by reabsorbing the factors in a redefinition of the spectral density as in lj21|l . 
This will additionally yield two boundary terms. 

In order to rewrite the other terms of l|58l) in terms of the sharp path, it is more convenient to rewrite relation 1|59|) 
in time domain 



y sm ,A(t') = / di'Vt 4 '- 4 ") [ 7 y S H,x(t") + (A£ - A£ c )6(t" - t )} . 



(62) 



Collecting all terms, we get the still exact result: 



P(\,t) 



E 



Tr R { p(t ) exp 



it) 



A£ + A )q 



.^(A£-A&) + (1- £ )A) V - 

h J r) 



iF 



dt'y S h,\(t') - 



(Ag - Age) 2 rTB 

2^ ^ R 
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iF(Ag - Ag e ) 
^7 



dt'y sh! x(t')Ll B (t' -t ) 



2h 



-A£ + A£ c -A 



(63) 



The influence phase ^py [y s h,A, £sh, i] is defined as in l|10(l . with the correlation functions L TB (r) and Mj TB (r) defined 
in terms of the new spectral density Jtb(w) given in Q21[) . Technically, 'J'™ [2/sh,Ai ^sh, i] can be rewritten in terms of 
the extended path y s h,\(t') by substituting everywhere y s h,\(t') for y s h,x(t') and extending the corresponding integrals 
to the whole real axis. This applies to both integrals in the real part, yielding a prefactor 1/2, and to the integrals 
on t' in the imaginary part. The average path x s h(t') was defined above Eq. (fT%|) . 

In order to go further, we now have to exploit the simplifications valid in the regime described at the end of 
Section [H] The long-time limit yields e <C 1 and the rare transitions limit A£ c -C A£. Furthermore, one can see 
that £™(0) diverges when the cutoff frequency uj c of an Ohmic spectral density J(lu) — ijLue^"^" tends to infinity, 
as in the strictly Ohmic case considered here. Then, in the configuration sum, the contributions of the configurations 
of {a} and {a'} for which the prefactor (Ag — A£ e ) 2 /2h r y 2 of L™(0) is minimal can be made as large as one will 
with respect to the contributions of other configurations. Therefore, combining with the rare transitions limit, we 
can restrict the configuration sum to configurations such that A£ = 0. This restricted sum is denoted with the boxed 
primed sum introduced in (|20[1 . Combining these properties, we find the much simpler expression 



P(X,t) 



]T' Tr R [p(t )e^ + ^} exp j ^[y sh ,x, x sh , t] 



iF 



dt'y sh , x (t') - 



ih\ 2 

2r] 



(64) 



This expression resembles already much more the one of a tight-binding generating function lj23|l . 

E. Extraction of the A-dependence 

The last task in order to demonstrate (|17|l is to extract explicitly the A-dependence out of the path y s h,x{t'), 
using H60f) . For the influence phase, we find 



$Fv[ysh,z s h,i] + A 



ft ft 



dt'y sh (t')N(t')-A x + dt'x sh {t')e 

V J t J to 



2rf 



N(t) - —N(t ] 
7 



with the auxiliary function N(t') — J. dt"N(t") given by 



N(t') 



dt" 



6(t" -t )-9(t" -t)-6(t" -t ) 



1-E 



Ll B (t'-t"). 



(65) 



(66) 



Using the explicit expression of the path y s h(t') given by (fTHf> and the property y s h{to) = — A£ = satisfied by the 
paths over which the configuration sum runs, one has 



dt'y sh {t')N(t') = L 



to 



x — > n x — ^ n 



In the rare transitions limit, it is possible to show 

N(t tI ) ~ 



2r]k B T 



ttT ~ ^ 



7 



when t tv equals any of the transition times tj , t'j, . Therefore one may rewrite 

t dt'y sh (t')N(t') = ^^ fdt'y sh {t'). 

t n Jt 



(67) 



(68) 



(69) 



Similarly, the term 



to 



En 
J'=l 



cr,e 



-7(<-* 3 ) 



E" 



(7/e 



' -7(*-</) 



is negligible in the rare transitions limit. Finally, for the driving contribution we obtain 



f , / , a & f i , , a AP 



t-to- 



1-e 



(70) 



(71) 



The last term in the square brackets may be neglected in the long-time limit. Putting these results in l|t)4|) completes 
the proof of the duality relation (|17|l . 



F. The tight-binding generating function at long time 

Here we shall prove that the solution of the generalized master equation Ij29(l behaves at long time as the expression 
given in l|3U|) . In order to discuss the long time behavior of the generating function, it is convenient to introduce its 
Laplace transform 



CPtb(X,s) 



dte- s (*- to )p TB (A,£). 



(72) 



The long-time behavior is related to the behavior of the Laplace transform at small s. Rewriting the generalized 
master equation for the Laplace transform turns the convolution into a product 



CPtb(X,s) 



1 



Ptb (A, t ) + / dre- sr K(X, t)£P tb (A, a) 



(73) 



This equation may easily be solved. Remembering that the initial populations are Pj TB (io) = <5/.o, yielding 
Ptb (A, to) = 1, we obtain 



£Ptb(A,s) 



dTe- ST K(X,T) 



(74) 



Let us compare this expression with the Laplace transform of the announced solution 130JI 



£Ptb(A,s) 



E r m ( 



e \mL _ 1 



(75) 



Remembering that K(X, r) = J2m=-oc eXmL Km(T), one sees that the two expressions are equal in the limit of small s 
if one defines 



/>OC 

r m = / dTK m (r) 
Jo 



(76) 



and provided the property To = — J2 m ^o ^ m ^ s satisfied. This property describes the conservation of the total 
population and can be verified on the explicit expressions for the transition rates T m . 



G. Explicit expressions for the transition rates 



We shall now discuss the explicit expression of the transition rates T m in a dissipative tight-binding model, char- 
acterized by the Hamiltonian ll'L'l) coupled to a bath of harmonic oscillators with spectral density Jtb(^)- Real-time 
path integrals techniques yield expressions for the transition rates in terms of a pair of tight-binding paths. Their 
structure is closely related to the one of the generating function (|23|l . in a way sketched in Section HV Fl As discussed 
m Section tTTTl one can classify the different contributions to a given rate T m with respect to the number N of tran- 
sitions occurring in the pair of tight-binding paths. These contributions form the A^th order rate denoted The 
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FIG. 4: Left: Graphical representation of the pairs of tight-binding paths corresponding to some contributions to the second- 
order (a,b, and c) and third-order (d and e) transition rates, with the corresponding dependence on the couplings A m . Right: The 
couplings associated with the transitions in the forward path (horizontal) and backward path (vertical). 



total rate is the sum of these contributions T m = X)5v=2 ^mP • The different contributions to the iVth order rates can 
themselves be split in terms of the numbers n, respectively n' , of transitions happening in the forward path (|18a|) . 



respectively backward path l|18b|l . This means = J2n=o J2n'=o *H 



,N)T 



(n,n ) 

= {(Ti, . . . , cr„; <t[, a' n ,}, which 
satisfy the constraints a j ~ J^j'—i a 'j' ~ m - The charges Uj, respectively a',,, characterize the transitions in the 



Each contribution to ft'™ can be parameterized by a set of charges a™'™ ' 



forward path l|18a|l . respectively backward path l)18b(l In the tight-binding model l)22|l. they can take any positive, 
respectively negative, integer value, representing a transition to the neighbor of cor resp onding order to the right, 
respectively to the left. For given to, n and n', these contributions may be written [2.ll9| 



with the prefactor 



n'-l 



/oo _ _ / f°° \ ~ / C°° \ 
dT Il / d ^ II / d P'j' exp{$ FV + *} 
I ...... I 'OO j-l \J0 / y = 1 \J0 / 



(77) 
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the influence phase 
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and the driving term 



fc=i j=i 
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(79) 



* = 



iFL 



(80) 



These last two expressions have to be rewritten in terms of the time variables pj = tj+i — tj, p'y = tj, +1 — t'-,, and 

T = 77 y^j'—i t'j> ~ n y^j— i, If = 0, respectively n' — 0, there are no integrals on r and p^, , respectively r and pj. 
The twice- integrated bath correlation function is given in 

The expression (|77l) is valid up to third-order N = n + n' = 2, 3 only. At higher order N > 3, reducible contributions 
have to be subtracted Qj| ■ 

The problem is to find all configurations of pairs of tight-binding paths involving N transitions in total, and satisfying 

the constraints X^=i °j = Xw'=i °j' = TO > or m °ther words to perform the configuration sum on am" ^ in Eq. lT77|l . A 
graphical representation can help us in this task. Let us start from a two-dimensional lattice parameterized by pairs of 
integers (I, I'). Let us represent by an horizontal arrow of length to to the right a transition in the forward tight-binding 
path to the neighbor of order to to the right. Such a transition is characterized by a charge a = +m. A transition to the 
left neighbor of order to, characterized by a charge a = —to, is represented by an arrow of length to to the left. Likewise, 
we represent transitions happening in the backward path, characterized by a charge a' = +to, respectively a' = —to, 
by vertical arrows of length to, pointing downwards, respectively upwards, in our diagram. Together, the arrows 
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|Ai| 2 


-Q(r) 


— r 



TABLE I: The different contributions to the rates Fm" \ at second order n + n = 2, involving only transitions to the 
nearest neighbors, characterized by charges Oj,oy = ±1, and mediated by the coupling A±i. Each contribution is specified 

by am'" ' = {ci, • • • , o n ; a[, . . . , a-' n ,}, given in the second column. The last three columns give the explicit expression of the 
prefactor, the influence phase, and the driving term for the corresponding contribution. 



corresponding to the transitions happening in the pair of tight-binding paths will draw a path in the diagram. With 
this representation, it is now easy to find the paths satisfying the constraints X)j=i a j = 2j'=i °j> ~ m - K one starts 
from a diagonal site of our diagram, say at position (0,0), the constraints select those paths which end up in the 
diagonal site (to, to). Some examples are given in Fig.|H Let us now look in detail to the second-order and third-order 
contributions. 

(2) 

In order to understand the different contributions to the second-order rates T m , let us first look at those contribu- 
tions which involve only transitions to the nearest neighbors. In this case, the charges aj and a'y may only take the 
values +1 and —1, denoting transitions to the nearest neighbor to the right, respectively to the left. These transitions 
are mediated by the couplings A±i. In other words, we consider the situation described by the Hamiltonian 122(1 
with A m = for |to| > 2. All contributing configurations are listed in Table [I] From there, it is clear that there is 

(2) (2) 

only one contribution to T\ , and only one to r_-[, which read 

= l^L dTe -(L 2 /h)Q(r)±i(FL/h)r^ (g^ 

(2) . . . (2) (2) (2) 

From the numerous contributions to Tg , we can explicitly check the relation T = — required in Sec- 

tion lTVFl 

Let us now allow all couplings A m in the Hamiltonian 1(22(1 to be nonzero and ask again the question: which pairs of 
tight-binding paths involving two transitions contribute to the transition rates? In the graphical representation shown 
in Fig.0|and described above, those paths are represented by a succession of two arrows connecting the initial diagonal 
site (0, 0) with any final diagonal site (m, m). It is clear that, if the first transition in any of the two paths reaches a 
neighbor site at distance to, the second transition must also happen between sites at the same distance to, in order to 
end up back in the diagonal. At second order, there will thus be no contributions to the rates combining two different 
couplings A m and A m > . One can convince oneself that all 6 possible combinations for a given to are similar to the ones 

listed in Tabled with ±1 replaced by ±m in am'™ ■ They are proportional to |A m | 2 and give contributions to T^ m 
and r o 2 ^. One easily sees that the general expression for the influence phase is $fv = —m 2 (L 2 /K)Q(t), and for the 
driving term ^ = ±\m{FL/K)T. Thus, the contributions to the second-order transition rates Y m , for any to ^ 0, can 
be described by the expression 1(34(1 . For m = 0, the relation To = — J2 m ^o^ m can ^ e explicitly verified. 

(3) 

We shall now investigate some contributions to the third-order rates . Let us first restrict ourselves to the case 
where, in the Hamiltonian ((22(1 . only the couplings to the nearest neighbors A±i and next-nearest neighbors A±2 
are nonzero. With the help of the graphical representation shown in Fig. 01 one finds 24 contributions in the a-a' 
description. They are listed in Tabled They all involve two transitions to a nearest neighbor and one to a next-nearest 
neighbor, in other words twice A±\ and once A±2- We shall group them under the notation Fm' ) [112] in order to 
distinguish them from other contributions which will arise when we shall switch on higher-order couplings A m , |to| > 2. 

Let us first look at [112]. Changing the integration variable variable r into — r and using the property Q(— r) = 

(121 (21) 

Q*(t), one can see that the second contribution T 2 is the complex conjugate of the first one T 2 ■ Put together, 



p(n,n' ) 
J- m 




(n.n'l 




A 


$PV 

(units of -^-) 


(units of i#) 


F (2,l) 
1 2 


{+l,+l;+2} 


All2 


Q(pi)-2Q(r+ipi)-2Q(r-ipi) 


2r 


r (l,2) 
1 2 


{+2;+l,+l} 


All2 


Q*(pi) - 2Q(r + ipi) - 2Q(r - ipi) 


2r 


r (2,l) 
1 1 




-l,+2;+l} 


— A*12 


-2Q( Pl ) + Q(t + ipi) - 2Q(r - ipi) 


T- fpi 


r (l,2) 
1 1 


{+2,-1; +1} 


— A 112 


-2Q( Pl ) - 2Q(t + ipi) + Q(r - ipi) 


r + fpi 


{+l;-l,+2} 


— All2 


-2Q*(p' 1 ) - 2Q(r + ipi) + Q(r - ipi) 


r + fpi 


r (3,0) 


{+i;+2,- 


1 x 
1 i 


— All2 


-2Q*(p'i) + Q(r + ipi) - 2Q(r - ipi) 


r - fpi 


{+ 


l,+l,-2 


0} 


— All2 


Q(p 1 )-2Q(p 1 +p 2 )-2Q(p 2 ) 


pi + 2p 2 




{+l,-2,+l 


0} 


— All2 


-2Q(p 1 ) + Q(p 1 +p 2 )-2Q(p 2 ) 


Pl - P2 




{- 


2.+1.+1 


0} 


— All2 


-2Q(p 1 )-2Q(p 1 +p 2 ) + Q(p 2 ) 


-2pi - p 2 




{- 


1,-1, +2 


0} 


All2 


Q(p 1 )-2Q(p 1 +p 2 )-2Q(p 2 ) 


-pi - 2p 2 




{- 


l,+2,-l 


0} 


All2 


-2Q(pi) + Q(pi + p 2 ) - 2Q(p 2 ) 


-Pl + P2 


p(0,3) 


{+2,-1,-1 


0} 


All2 


-2Q(p 1 )-2Q(p 1 +p 2 ) + Q(p 2 ) 


+2pi + p 2 


{0 




-2} 


— A*12 


Q*(pi)-2Q*(pi+p 2 )-2Q*(p 2 ) 


-pi - 2p 2 




{0 


+l,-2,+l} 


— A 112 


-2Q*(pi) + Q*(pi+p 2 )-2Q*(p 2 ) 


-Pl +P2 




{0 


-2,+l,+l} 


— All2 


-2Q*(pi)-2Q*(pi+p 2 ) + Q*(p 2 ) 


2pi + p 2 




{0 


-1,-1, +2} 


All2 


Q*(pi)-2Q*(pi+p 2 )-2Q*(p 2 ) 


pi + 2p 2 




{0 


-l,+2, 


-1} 


All2 


-2Q*(pi) + Q*(pi + p 2 ) - 2Q*(p 2 ) 


Pi - P2 


r (2,l) 
1 -1 


{0 


+2,-1, 


-1} 


All2 


-2Q*(p'i) - 2Q*(pi + p 2 ) + Q*(p 2 ) 


-2pi - p 2 


{+1,-2;- 


1} 


All2 


-2Q( Pl ) + Q(r + ipi) - 2Q(r - ipi) 


-T+fpi 


r (l,2) 
1 -1 


{- 


-2,+l;- 


1} 


All2 


-2Q( Pl ) - 2Q(r + ipx) + Q(r - ipi) 


-T- fpi 


{- 


-i;+i,- 


2} 


All2 


-2Q*(pi) - 2Q(t + ipi) + Q(r - ipi) 


- fpi 




{- 


-i;-2,+l} 


All2 


-2Q*(pi) + Q(r + fpi) - 2Q(r - ipi) 


-T + fpi 


r (2,l) 
1 -2 


{- 


-1,-1;- 


2} 


— A*12 


Q(pi)-2Q(r+ipi)-2Q(r-ipi) 


-2r 


r (l,2) 
1 -2 


{- 


-2;-l,- 


1} 


— All2 


Q*(pi)-2Q(r+ip'i)-2Q(r-ipi) 


-2r 



TABLE II: The different contributions to the rates Tm' n \ at third order n + n' — 3, involving only transitions to the nearest 
neighbors and next-nearest neighbors, characterized by charges Oj, oy = ±1, ±2, and mediated by the couplings A±i and A±2. 

Each contribution is specified by a^'" ' = {ui, . . . , a n \ a[, . . . , o' n ,}, given in the second column. The third column gives the 
explicit expression of the prefactor, in terms of A112 = — iA? A^/h 3 , for the corresponding contribution. The last two columns 
give the influence phase and the driving term. 



they can be rewritten 



rf >[ii2] = Aim |a?a* j°° dr J™ d P e^^- 2 ^Hp)-2Q(T-y)]+^T 



Likewise, we obtain 



r (_ 3 )[ll2] = ^lm^AfA 2 J°° dr dpc^^-^+^-^-y^^Y 



(82) 



(83) 



Following the same procedure, one can show that the contributions to if' 1 - 1 are the complex conjugate of the ones 
to r^' 2 \ yielding 

if ) [112] = -Alm{A 2 A* J°° dT^° C dp[c*[- 2 «-'')- 2 ^+^)+^-^]+ i *^+i'') 

+ e #[-2Q(-p)+Q(r+ip)-2Q(r-Ip)]+i£i(r-|rtj j ^ 



Both parts of the expression can be partially merged by substituting r — |p for r in the first one, and t + |p for r 



in the second one. This yields 



rj 3) [ii2] 



h 3 



1ml A? A* 



dr 



dpe 



: Q(-p)+i J 



e * 



-2Q{r+p)+Q(T+2p)} 



e » 



\-2Q{T-p)+Q(T-2p)} 



(85) 



A similar expression is obtained for 



rL 3 ![ii2] 



dr 



dpe 



-Q(-p)-U 



i-[-2Q(r+p)+Q(T+2p)] + e ^[_2Q(T-p)+Q(r-2p)] 



Upon using the modulus and phase of the couplings A r , 
common expression (136(1 . 



I A. 



all these rates may be rewritten in the 



The rate Tq 3 ^ [112] here again satisfies T^' [112] = — J2 m ^o ^m* [112]. This relation can be explicitly checked with the 
expressions given in Table ITT1 One has to take into account the various chronological orderings of the transition times 
in the forward and backward paths, which lead to 48 configurations. By proper substitutions in the time integrals, 
one can then check that they compensate each other two by two. 

Let us now allow the couplings A±3 to be nonzero. Still using the graphical representation shown in Fig. one 
finds 48 additional contributions to the third-order rates. They all involve one transition of each of the three different 
kinds, therefore we group them under the notation Fm [123]. Applying the same procedure as for the third-order 
contributions involving only A±i and A±2 7 we arrive at the expression l|38[l. 

The generalization to higher-range couplings A m , |m| > 3, and larger number of transitions now becomes clear. If we 
would allow the couplings A±4 to be nonzero as well, we would get additional third-order contributions involving A±4 
with A±3 and A±i, as well as contributions involving A±4 and twice A±a- hi order to write down the explicit 
expressions for these contributions, one has to follow a procedure similar to what we did here above. Going to 
fourth or higher order TV, the procedure would also remain the same. One can see that the number of contributions 
increases factorially with the number of transitions. Starting from fourth order, one has also to subtract the reducible 
contributions in order to evaluate the transition rates. 



,(3) 



(3)r 



V. CONCLUSION 



We have developed a method yielding the duality relations (|17J) and ((26(1 between the long-time dynamics in a 
tilted ratchet potential in the presence of dissipation, and the long-time dynamics in a driven dissipative tight-binding 
model. Detailed proofs have been presented. The formalism has been applied to the evaluation of the current in 
quantum ratchet systems, yielding an expression ((33(1 in terms of the transition rates in the tight-binding system. 
Consequently, the evaluation of these rates has been discussed in Section IlIII In particular, the obtained expressions 
show the explicit dependence of the ratchet current on the parameters of the ratchet potential. 

This approach allows to investigate quantum ratchet systems in the weak dissipation limit, which is beyond the 
validity range of many of the other theoretical approaches. With respect to the perturbative approach of Ref. Q, the 
method reported here has the advantage that the nonlinear regime of a large driving force, eventually leading into the 
classical regime, can be reached. On the other hand, the continuous system considered is essentially different from 
the tight-binding molecular wire investigated in Ref. 6], and may thus apply to different experimental situations. 
Weak dissipation is even a favorable situation in our approach. Indeed, the duality relation links a situation of weak 
dissipation in the original model with strong dissipation in the dual tight-binding model, and vice versa. Therefore, 
one is brought to evaluate the tight-binding transition rates in the limit of strong dissipation, where it suffices to 
consider the lowest orders in the tunneling amplitude. 
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